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Abstract 

In this article we propose a novel face recognition method based on Principal 
Component Analysis (PCA) and Log-Gabor filters. The main advantages of the 
proposed method are its simple implementation, training, and very high recognition 
accuracy. For recognition experiments we used 5151 face images of 1311 persons 
from different sets of the FERET and AR databases that allow to analyze how 
recognition accuracy is affected by the change of facial expressions, illumination, and 
aging. Recognition experiments with the FERET database (containing photographs 
of 1196 persons) showed that our method can achieve maximal 97-98% first one 
recognition rate and 0.3-0.4% Equal Error Rate. The experiments also showed that 
the accuracy of our method is less affected by eye location errors and used image 
normalization method than of traditional PCA -based recognition method. 

Key words: Face recognition, Principal Component Analysis, Log-Gabor filters, 
FERET database 



1 Introduction 



Principal Component Analysis (PCA) or Karhune n Loeve Transform (KIT) 
- based face recognition method was proposed in ((Turk and Pentlandl . Il99lh 
and became very popular because of its relatively simple implementation 
and high recognition accuracy. During past fifteen years face recognition was 
a field of active research, and many other statistical methods (related to 
PCA) were investigated and proposed for face recognition: Line ar Discriminant 
Analysis (LDA), Independent Component Analysis (ICA) ( Bartlett et all 
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20021) . Kernel PCA, Dual PCA (|Moghaddairl 120021 ). Because KLT is data 
dependent and is not very fast, other transform s were also used for fac e 
recognition: Discrete Cosine Tr ansform (DCT) ( F afed and Levind . 2001), 
Fast Fourier Trans form fFFT) fapi^ M d Rj^ bj , Discrete Wavelet 
Transform fDWTl (iFeng et all EfflOl) . Wave let Packet Decomposition (WPD) 
f Garcia et all 120001 )! IWiskott et al. (1997) proposed Elastic Bunch Graph 
Matching (EBGM) and Gabor wavelets -based face recognition m ethod that 
achieved very high recognition accuracy. Escobar and Solar (|2002l, used EBGM 
-based recognition of faces in Log-Polar coordinates. C omprehensive over view 



of various face recognition methods could be found in ( Zhao et al. . 2000[ ). As 



it was shown by numerous experiments, face recognition accura cy can be in- 



creas ed by combinin g several methods, for example, D CT+PCA (IRamasubramanian and Venkatesh . 
200lh DWT+PCA (jFeng et al!l2000h . WPD+PCA (|Perlibakasl . l2004l ). PCA+LDA 



(jZhao et"aL 1998f ). or by using various image pre-processing methods. Recent 
results also showed that using Gabor or Log-Gabor features instead of tra- 
ditional greyscale features and by combining these features with well known 
recognition methods like PCA, ICA, LDA or SVM it is possible to achieve 
very high recognition acuracy. Now we will overview various combined face 
recognition methods that use Gabor or Log-Gabor features and that are re- 
lated with a method that we propose in this publication. Although many 
researchers used the same Gabor filters and well known feature compression 
and classification methods, all proposed methods differ from each other by 
feature selection techniques, parameters of filters, used classification method 
and its parameters, distance measure, and image normalization method. Af- 
ter filtering face image with Gabor filters of multiple scales (usually 4-5) and 
orientations (usually 6-8) we get very large number of features (24-40 images 
of the same dimensions as initial image). Perhaps one of the most important 
questions is how to reduce the number of these features for further process- 
ing. The most popular method is to extract Gabor features at a small number 
(usually less than 100) of face poin ts around face featur es (like eyes, lips, nose) 



that were detected using EBGM ((Wiskott et all 119971 ) or similar method. At 



each d etected point are ex tracted Gabor features from all scales and orienta- 



tions. ILvons et al. ( 2000l) combined EBGM Gab or features and LDA-based 
recognition method. Smeraldi and Bigun ( 20021 , detected face features us- 
ing saccadic search with a set of Log-Gabor filters that were arranged to 
concentric circles (retinas). At detected points were extracted Log-Gabor fea- 
tures and passed t o the SVM -based classi fier. EBGM -based feature selection 
was also used by IWang and Tangl (2003) for Bayesian PCA -based recogni- 
tion. Because EBGM -based methods require training with manually labelled 
faces, other researchers use feature extraction methods that do not require 
such training. One possible approach is to combine Gabor magnitude images 
from all scales and orientations to a single feature vector (image) and use this 
vector for recognition. Such method was used by Zhang et al. ( 2004 ) for Ga- 
bor+AdaBoost face recognition. Because such combined feature vectors may 
be too large for further processing, we can reduce the number of features by 
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using smaller initial images. Liu and Wechslei ( 2002h decided to downsample 
feature images of each scale and resolution a nd then combine these fe atures to 
a single vec tor for Gabor+Enhanced LDA ( Liu and Wechsler . 20021 and Ga- 
bor+ICA ( Liu and Wechsler . 2003j ) -based recognition using L\, Euclidean, 
and cosine -b ased distance measures . In order to reduce the number of Ga- 
bor features, iKepenekci et al. (2002) used sliding window based search at all 
scales and orientations. In each window were extracted features with maximal 
magnitudes, stored their locations, and then both magnitudes and locations 
were used for comparison . For recognition was used very similar distance mea- 
sure that was used by ([Wiskott et all Il997l ) with additional constraints to 



feature locations. 



In this article we propose to find the locations of Log-Gabor (jFieldi Il987l) 
features with maximal magnitudes at single scale and multiple orientations 
using sliding window -based search and then use the same feature locations 
for all other scales. For further feature compression we used Principal Com- 
ponent Analysis (PCA) because its simple implementation, fast training and 
because using PCA with "whitened" angle -based distance measure it is possi- 
ble to achieve similar r ecognition a ccuracy like using EBGM and LDA -based 



recog nition methods ijCSlA 120031 ). We tested our method using 5151 face 



images of 1311 persons from the FERET and AR databases and the results 
showed that the proposed recognition method can achieve higher recognition 
accuracy than many other existing methods. The results of experiments also 
showed that PCA using Log-Gabor features is less sensitive to face detec- 
tion errors and used image normalization method than PCA using greyscale 
features. 



2 Feature extraction using Log-Gabor filters 



The Log-Gabor filters were proposed by iFieldl ()1987f ) for coding of natural 
images. The experiments showed, that these filters are consistent with the 
measurements of the mammal ian visual syst em and are more suitable for cod- 
ing of natural images than iGabor (1946) filters. The Log-Gabor filter in 
frequency domain can be constructed in terms of two components, namely the 
radial filter component G(f) and the angular filter component G(8). In polar 
coordinates the filter transfer function co ul d be wr i tten i n the following form 



( Fieldlll987h . ijBienin and du Bufl . ll994h . (|Kovesil . Il996) 



G(f,9)=G(f)-G(9) = ex P 



(log(///o)) 2 * 
2(log(k/f ))\ 



■ exp 



2ai 



here /o is the centre frequency of the filter, k determines the bandwidth of 
the filter, o is the orientation angle of the filter, and a e = A 9/sq where sq - 
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scaling factor, A 9 - orientation spacing between filters. For face recognition we 
generated multiple Log-Gabor filters of different scales and orientations using 
the following parameters: f = 1/A, A = A • s^ s X \ kj f = <Jf, n s = 1, N s ; 
6 = ix{n - 1)/N a , A 9 = n/N , n = 1, N Q ; A = 5, s x = 1.6, a f = 0.75, 
N s = 4, sg = 1.5, N Q = 6, here Ao is the wavelength of the smallest scale filter, 
s\ is the scaling factor between successive filter scales, N s is the number of 
scales, N a is the number of orientati ons. Most of t hese parameters were chosen 
following the recommendations of ( Kovesil . 2003| ). 



Using Eq. 1 we calculate two-dimensional Log-Gabor filter G nojTls in Fourier 
space of a chosen filter scale and orientation. The size of the filter array G no ^ ng is 
the same as the size of the two-dimensional image / that we wish to filter. Then 
we perform filtering (convolution in Fourier space), magnitude calculation and 
masking using the following equation: 

V no , ns = abs(IFFT2(G no , ns . * FFT2(I))). * mask, (2) 

here ".*" - array (not matrix) multiplication, / - normalized (cropped, masked) 
face image, G no , ns - Log-Gabor filter of desired orientation and scale in Fourier 
space, FFT2 - two-dimensional Fast Fourier Transform, IFFT2 - inverse 
FFT2, mask - binary mask for masking magnitude image (the same as is 
used for masking greyscale face image I in order to leave only the internal 
part of the face), V notTls - masked Log-Gabor magnitude image. 

After image filtering with multiple Log-Gabor filters (iV s scales and N Q orien- 
tations) we get very large number of Log-Gabor features (magnitude values 
in all N s ■ N a magnitude images). In order to reduce the number of features 
and achieve partial face recognition invariance with respect to different facial 
expressions and minor face detection errors, we use sliding window algorithm 
that is illustrated in Fig. 1. Rectangular window of a chosen size (e.g., 8x8 
pixels) is slided over the magnitude image V noj i using some sliding step (e.g., 
6 pixels, overlapping of windows is 8-6=2 pixels). In each window we find one 
maximal magnitude value and remember the location (coordinates in image 
V A „ o> i) of this value. If several equal values are found, we use the one that is 
closer to the centre of the window. If magnitude image is masked, we perform 
search only in an unmasked image part. 

We apply sliding window algorithm only for the first scale (n s = 1) (obtained 
using filter with the smallest chosen wavelength) of each orientation and find 
the locations of highest magnitudes. Then these locations are used for extract- 
ing magnitudes from other scales, corresponding to the analysed orientation, 
as it is illustrated in Fig. 1. In this figure feature locations are marked with 
black points. Using scales n s = 1 we decide at what locations (coordinates) 
we will extract the features and then use the same coordinates for all mag- 
nitude images, corresponding to the same processed orientation n Q (the same 
orientations n Q - the same locations, different orientations - different loca- 
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tions). All extracted Log-Gabor features (magnitude values) are stored in a 
one-dimensional vector X and used as input for Principal Component Analysis 
-based face recognition method. 




Fig. 1. Selection of Log-Gabor magnitude features using sliding window algorithm 

Example Log-Gabor magnitude images are presented in Fig. 2 (images are 
inverted, dark points mean high magnitude values). In this example we used 
filters of N = 6 orientations and N s = 4 scales (N a ■ N s = 24 filters), and 
filtered a normalized (derotated, masked) facial image with these filters. Cal- 
culated Log-Gabor magnitude images were also masked. Left-most binary im- 
ages in Fig. 2 show the locations (black points) of Log-Gabor features that 
were found using sliding window algorithm. It must be noted, that for the 
selected image size the Log-Gabor filters (of different sizes and orientations) 
can be calculated only once and stored. When we perform face recognition, 
the Log-Gabor features (found using sliding window) for each image from the 
database of faces are also calculated only once and stored. 



3 Face recognition using Principal Component Analysis of Log- 
Gabor features 

In this section we will describe Karhunen-Loeve transform (KLT) -based face 
recognition method, that is often called Principal Component Analysis (PCA). 
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Fig. 2. Image magnitudes after using Log-Gabor niters, and the locations of features 
at different orientations that were found using sliding window algorithm 

We will prese nt on l y the main formulas of this method, which details could 
be found in (|QrofiL Il994h . 



Let Xj be iV-element one-dimensional image-column (vector) and suppose 
that we have r such images (j = l,...,r). In traditional PCA -based face 
recognition method, these images contain grey values of the two-dimensional 
facial photographs. In our case these one- dimensional images Xj (data vec- 
tors) contain Log-Gabor features. We calculate the mean vector, centred data 

r r 

vectors and covariance matrix: m = - J2 Xj, dj = Xj — m, C = - J2 djd 



here X = (xi, X2, ...,xn) t , rn = (mi, m 2 , mjv) T , d = (di, d^ 



d 



Nj 



3=1 



r 

3 ' 



In order to perform KLT, it is necessary to find eigenvectors and eigenvalues 
Afc of the covariance matrix C ( Cuk = A^u^ ). Because the dimensionality (N 2 ) 
of the matrix C is usually large even for small images, and computation of the 
eigenvectors using traditional methods is complicated, di mensionality of ma- 
trix C is reduced using the decomposition described in ((Kirbv and Sirovich . 
1990) (if the number of training images is smaller than the length of the vector 
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X). Found eigenvectors u = (ui, u 2 , ...,un) t are normed and sorted in decreas- 
ing order according to the corresponding eigenvalues. Then these vectors are 
transposed and arranged to form the row-vectors of the transformation ma- 
trix T. Now an y data X can be projected into the eigenspace and "whitened" 
( Bishopl . 1995| ) using the following formula: 



Y = A- 1/2 T(X -m), (3) 
here X = (x u x 2 , .-,x N ) T , Y = (y u y 2 , y r , 0, 0) T , 



A-V2 = diagiJl/X!, Jl/)C 2 , sJl/K). 

For projection we can use not all found eigenvectors, but only a few of them, 
corresponding to the largest eigenvalues. We can manua lly select the desir ed 
number of eigenvectors or use the method described in ( Swets et al. . 19981 ). 



For each facial image (that we wish to use for recognition) we find Log-Gabor 
features, projected these features into the eigenspace and calculate eigenfea- 
ture vector Z = (zi, z 2 , z n ) T = (yi, y 2 , y n ) T , here n is the number of 
features. Recognition of unknown face is performed by calculating its feature 
vector Z new and comparing it with the feature vectors of known faces. For com- 
parison we calculate the distances Ei(Z new , Zi) between unknown face and each 
known face and say that the face with feature vector Z new belongs to a person 
s = arg min . For rejection of unknown faces a threshold r is chosen and 

i 

it is said that the face with projection Z new is unknown if e s > r. For recog- 
nition we used cosine-based distance measure Ei(Z new , Zi) = —cos(Z new ,Zi), 
because using this distan ce measure we can achieve higher recognition ac- 
curacy ( Perlibakasl . 2004 ) than using the Euclidean or Manhattan distance 



measures. 



4 Normalization of face images 



For recognition experiments we used two image normalization methods. One 
method uses manually selected centres of eyes and the tip of chin (3-point 
normalization method), and another method for normalization uses only the 
centres of eyes (2-point normalization method). Image normalization proce- 
dure of 3-point method is presented in Fig. 3. The last image (Fig. 3 (e)) 
also presents the result of 2-point normalization. For illustration we used an 
image from our personal archive. 

Now we will describe our implementation of 3-point normalization method. 
Initial images were denoised (using Gaussian filter with a = 0.5 and window 
size 5x5), derotated (in order to make the line connecting eye centres hori- 
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(a) (b) (c) (d) 

Fig. 3. Image normalization (a-d - 3-point method, e - 2-point method): a) ini- 
tial image with selected eyes and chin b) denoised, derotated image, and cropping 
schema c) cropped, resized and masked image d) normalized image after histogram 
equalization 

zontal), cropped, resized (to the size of 128x128 pixels), masked. For rotation 
and resizing we used bicubic interpolation. For masking we used an ellipse 
with central point (64.5,45.5), horizontal axis of 120 pixels, and vertical axis 
of 160 pixels. Then for unmasked part of the image we performed histogram 
equalization (256 levels). When the image is masked, are left 12646 unmasked 
pixels of 16384 (128x128). 



For initial comparison of PCA and Log-Gabor PCA methods we used 3-point 
normalization method in order to perform experiments with faces that are not 
overcropped and also contain no scene's background information. Because the 
tip of chin may be hard to locate, most recognition methods for normalization 
use only the centres of eyes. So we used 3-point normalization only for initial 
comparison of PCA and Log-Gabor PCA methods, and for the rest of ex- 
periments we used 2-point normalization method (for normalization are used 
only the centres of eyes). Because there is no agreement how images should 
be normalized for face recognition experi ments , we im plemented 2-point nor- 
malization that is very similar to the (jCSll [2003) normalization method. 



Simil ar method was also used by some participants of the FERET f NISTl . 



2001) tests. At first images are derotated in order to make the line connecting 



eye centres horizontal. Then images are resized in order to make the distance 
between eyes equal to 70 pixels and cropped to the size of 130x150 pixels. 
During cropping the centres of eyes are vertically positioned on y=45 line (the 
centre of coordinates (0,0) is in the left top corner). Then the image is masked 
using an ellipse with its central point (65.5,50.5), horizontal axis of 128 pixels, 
and vertical axis of 236 pixels. After masking are left 17237 unmasked pixels. 
For an unmasked image part is performed histogr am equalizat ion. The main 
differences between our 2-point normalization and ( CSUl . 2003h normalization 
are as follows: initial images we filtered using Gaussian filter (CSU used no 
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filtering), for image rotation and resizing we used bicubic interpolation (CSU 
used bilinear interpolation), cropped 130x150 images we resized to 128x128 
pixels and then masked with resized (to 128x128) binary mask. After such 
masking were left 14454 unmasked pixels. Then for an unmasked image part 
we performed histogram equalization. The result of this normalization is pre- 
sented in Fig. 3 (e). 

When we performed face recognition using Log-Gabor PCA method, masked 
face images were filtered with Log-Gabor filters (24 filters of 6 orientations and 
4 scales) and calculated magnitude images. We masked these Log-Gabor mag- 
nitude images using the same masks that were used for image normalization 
and performed sliding window search of Log-Gabor features. Search window 
size is 8x8 pixels, sliding step is 6 pixels (the same in horizontal and vertical 
directions), and window overlap is 8-6=2 pixels. After using sliding window 
algorithm with masked magnitude images, we select 9240 magnitude values 
(Log-Gabor features) for 3-point normalization method and 10008 magnitude 
values for 2-point normalization method, that are located in the unmasked 
parts of magnitude images. This is the total number of values in all 6 orien- 
tations and 4 scales. Also we can use an unmasked magnitude images (initial 
greyscale images are always masked) , perform sliding window search in a whole 
magnitude image and select in total 10584 features for both 3-point and 2- 
point methods (the size of images is the same). In all the experiments we used 
the same normalized image patterns and the same implementation of PCA. 
The distances between "whitened" feature vectors were measured using cosine 
-based distance measure. 



5 Used recognition performance measures 



For comparison of face recognition methods we used Cumulative Match Char- 
acteristic (CMC) a nd Receiver Op erating Characteristic (ROC) - based mea- 
sures described in ((Brombal . 2003^ 1: the area above Cumulative Match Charac- 
teristic (CMCA) (smaller CMCA means better overall recognition accuracy); 
how many images (in percents) must be extracted from the database in order 
to achieve some cumulative recognition rate (e.g., not smaller than 95-100%) 
(smaller values mean that we need to extract fewer images in order to achieve 
some cumulative recognition rate); Equal Error Rate (EER) and the area be- 
low Receiver Operating Characteristic (ROCA) (smaller values mean better 
results); first one recognition rate (First 1) that is achieved if only the first one 
(most similar) image from the database is extracted (larger values mean better 
result). Percent (rank) of images that we need to extract from the database in 
order to achieve 100% cumulative recognition rate in the future we will denote 
as CumlOO. Graphical representation of the used characteristics is shown in 
Fig. 4- 5. 
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Fig. 4. Cumulative Match Characteristic Fig. 5. Receiver Operating Characteristic 
(CMC) (ROC) 



6 Experiments and results 



For r ecognition experiments we used the FERET database (Ph illips et al.. a . 
199£jj) containing greyscale photographs of 1196 persons. This database was 
collected in 1993-1996 at George Mason University during the FERET (FacE 
REcognition Technology) program. As far as we know, this is one of the 
largest databases of face photographs (of different persons') in the worl d that 



is publicly available for face recognition research purp oses (lQrofll2005l ) This 



database is widel y used for evaluati ng identification ((Phillips et ah. al 1 



and verification ( Rizvi et. all 119981 ) performance of face recognition methods. 



For recognition experiments we used 3541 facial images from this database. 
The size of each image is 256x384 pixels, for each image this database contains 
manually selected eye coordinates. For training we used 1196 greyscale images 
from the fa set of this database. The same 1196 fa images were used as a 
gallery (known persons), and images from other sets (1195 fb images, 722 
dupl images, 234 dup2 images, 194 fc images) were used as probes (unknown 
persons that we wish to recognize), fb set contain face images with different 
facial expressions, dupl and dup2 sets contain images that were taken after 
some time interval (up to 1.5 years) from fa images, and fc set contains images 
with different image input conditions (camera position and illumination). 

At first we performed recognition experiments with fa and fb sets using dif- 
ferent number (100-1000) of PCA features (different number of used eigen- 
vectors, corresponding to the largest eigenvalues) and compared the proposed 
masked Log-Gabor PCA with traditional PCA (using cosine -based distance 
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measure between "whitened" feature vectors). Images were normalized using 
3-point normalization method. The results are presented in Table 1. The 
results showed, that first one recognition rate of masked Log-Gabor PCA 
(89.29-98.41%) is always higher than of traditional PCA (80.42-88.03%), and 
EER values of masked Log-Gabor PCA (0.33-1.59 %) are always lower than 
EER values of traditional PCA (1.92-4.26 %). Other characteristics (CMCA, 
ROCA, cumulative recognition) of masked Log-Gabor PCA are also better 
than of traditional PCA. The results showed, that masked Log-Gabor PCA 
achieves larger first one recognition accuracy when we use larger number of 
features (e.g., 100 PCA features - 89.29% recognition accuracy, 1000 PCA 
features - 98.41% recognition accuracy). It must be noted, that masked Log- 
Gabor PCA uses shorter vectors (9240 Log-Gabor features) than traditional 
PCA (12646 greyscale features) for PCA training. Also we investigated an- 
other version of Log-Gabor PCA when Log-Gabor magnitude images are not 
masked (only magnitudes, initial images are always masked). In this case the 
sliding window search selects 10584 Log-Gabor features. The results showed, 
that in some cases unmasked Log-Gabor PCA can achieve higher first one 
recognition accuracy, but because these differences are not very large (<0.2% 
with >200 features) we prefer to use masked version of Log-Gabor PCA. 



Also we compared our face recognition results with the results of other re- 
searchers. For comparison we used the best results of the FERET program 
participants t hat took official FERET 1996-1997 tests (|Phillips et al.U2000h . 
f NISI! 2001 ). Also we present the best results of some other researchers 
who tested their recognition methods using all images (not subsets) from the 
FERET fa (images of 1196 persons) and fb (images of 1195 persons) sets 
that contain faces with different facial expressions. Most part of the compared 
methods for face normalization used manually located coordinates of eye cen- 
tres. These coordinates were marked by the creators of the FERET database 
and are distributed with this database. The results are summarized in Table 
2. Fully automatical methods are denoted by "auto", our different normal- 
ization methods are denoted by "3 pt." (3-point normalization) and "2 pt." 
(2-point normalization). In Figures 6 - 7 we also present CMC and ROC 
characteristics of our Log-Gabor PCA method (900 PCA features, 2-point 
normalization) . 



Now we will briefly describe face recognition methods of other researchers and 
will compare achieved results. MIT 1996 (Massachusetts I nstitute of Technol- 
ogy, MIT Media Laboratory) method was developed by ([Moghaddam et aT 



1996). For recognition they used dual (intrapersonal and extrapersonal) PCA 



and Bayesian MAP (maximum a posteriori) similarity measure. For learning 
were used image pairs of the same and different persons. UMD 1996, UMD 
1997 (University of Maryland) face recognition methods are based on PCA and 
LDA (Linear Discrimin ant Analysis) and were developed by ((Etemad and Chellapa , 
1997t[zhao et al.l . ll998| ). For training were used several images per person, for 
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Table 1 

Comparison of 3 face recognition methods: 1) Masked Log-Gabor PCA (MskLG), 
2) UnMasked Log-Gabor PCA (UnMskLG), and 3) Traditional PCA (Trad). 



Method Feat. Rank (%) in order to achieve desired CMCA First 1 EER, ROCA 
num. cumulative recognition accuracy rec, % % 







95 


96 


97 


98 


99 


100 










MskLG 


100 


0.33 


0.42 


0.50 


0.84 


2.34 


18.14 


17.67 


89.29 


1.59 


9.49 


UnMskLG 


100 


0.42 


0.50 


0.75 


1.17 


2.17 


11.87 


17.20 


89.87 


1.42 


9.03 


Trad 


100 


0.84 


1.17 


1.59 


2.93 


6.44 


36.87 


37.27 


83.85 


2.26 


28.77 


MskLG 


200 


0.17 


0.25 


0.33 


0.42 


0.84 


6.69 


11.34 


93.56 


0.75 


2.48 


UnMskLG 


200 


0.17 


0.17 


0.25 


0.42 


0.75 


6.94 


11.44 


93.56 


0.75 


2.61 


Trad 


200 


0.59 


0.84 


1.17 


2.34 


7.19 


84.20 


40.21 


86.44 


2.01 


31.59 


MskLG 


300 


0.08 


0.17 


0.17 


0.25 


0.42 


7.27 


10.36 


95.31 


0.59 


1.65 


UnMskLG 


300 


0.08 


0.17 


0.17 


0.25 


0.59 


8.61 


10.41 


95.40 


0.50 


1.73 


Trad 


300 


0.59 


0.84 


1.34 


2.84 


5.52 


93.65 


39.36 


87.87 


1.92 


29.91 


MskLG 


400 


0.08 


0.08 


0.17 


0.17 


0.33 


3.76 


9.46 


96.99 


0.42 


0.93 


UnMskLG 


400 


0.08 


0.08 


0.08 


0.17 


0.33 


6.44 


9.58 


97.15 


0.42 


1.04 


Trad 


400 


0.59 


0.92 


1.67 


3.09 


12.12 


92.73 


52.94 


88.03 


2.09 


42.79 


MskLG 


500 


0.08 


0.08 


0.08 


0.17 


0.33 


2.84 


9.39 


97.24 


0.42 


0.69 


UnMskLG 


500 


0.08 


0.08 


0.08 


0.17 


0.33 


3.68 


9.26 


97.15 


0.42 


0.63 


Trad 


500 


0.92 


1.25 


2.01 


4.10 


16.64 


74.50 


58.73 


87.53 


2.34 


48.26 


MskLG 


600 


0.08 


0.08 


0.08 


0.17 


0.25 


2.26 


9.07 


97.74 


0.33 


0.48 


UnMskLG 


600 


0.08 


0.08 


0.08 


0.17 


0.25 


1.42 


8.91 


97.74 


0.33 


0.39 


Trad 


600 


1.34 


1.67 


2.59 


5.94 


13.55 


82.86 


67.05 


86.03 


2.59 


56.21 


MskLG 


700 


0.08 


0.08 


0.08 


0.17 


0.17 


2.76 


8.89 


97.91 


0.42 


0.36 


UnMskLG 


700 


0.08 


0.08 


0.08 


0.08 


0.17 


2.42 


8.86 


98.08 


0.33 


0.32 


Trad 


700 


1.42 


2.26 


4.35 


8.86 


22.41 


95.07 


88.53 


85.02 


3.18 


77.03 


MskLG 


800 


0.08 


0.08 


0.08 


0.17 


0.17 


3.51 


8.98 


97.99 


0.33 


0.38 


UnMskLG 


800 


0.08 


0.08 


0.08 


0.08 


0.17 


3.18 


9.02 


98.08 


0.33 


0.38 


Trad 


800 


2.09 


3.34 


5.35 


10.28 


28.51 


93.14 


109.44 


83.09 


3.51 


97.12 


MskLG 


900 


0.08 


0.08 


0.08 


0.08 


0.25 


1.84 


8.92 


98.16 


0.33 


0.33 


UnMskLG 


900 


0.08 


0.08 


0.08 


0.17 


0.17 


1.76 


8.85 


98.08 


0.33 


0.28 


Trad 


900 


2.76 


4.35 


7.02 


12.63 


38.13 


92.89 


123.40 


82.43 


3.93 


110.99 


MskLG 


1000 


0.08 


0.08 


0.08 


0.08 


0.17 


6.35 


9.44 


98.41 


0.33 


0.63 


UnMskLG 


1000 


0.08 


0.08 


0.08 


0.08 


0.17 


2.51 


9.01 


98.24 


0.33 


0.35 


Trad 


1000 


4.01 


5.52 


9.28 


16.22 


45.48 


90.80 


143.82 


80.42 


4.26 


130.80 
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Table 2. Recognition of expression -variant faces from the FERET database (gallery contains 1196 fa images, and probe contains 1195 
fb images). 



Method and its 




Rank (%) in 


order to achieve 




CMCA, 


First 1 


EER, 


ROCA, 


authors 


desired cumulative 


recognition, (0, 100%] 


[0, 10 4 1 

L " J 


recognition, 


[0, 100%] 


[0, 10 4 ] 

L " J 




95 


96 


97 


98 


99 


100 




[0, 100%] 






MTT 1996 fMoffhaddam et ah. 19961 


0.17 


0.25 


0.33 


1.09 


23.33 


99.83 


84.14 


94.81 


4.77 


203.25 


MIT 1996 auto fPhillios et ah. b. 1998) 
















-88.00 






UMD 1996 CEtemad and Chellaoa. 1997) 
















-83.50 


-7.00 




UMD 1997 fZhao et ah. 1998) 


0.08 


0.08 


0.17 


0.33 


0.84 


75.92 


18.91 


96.23 


1.09 


14.37 


USC 1997 (Okada et ah. 1998) 


0.17 


0.17 


0.25 


0.33 


3.09 


50.67 


27.94 


94.98 


2.51 


57.52 


USC 1997 auto fOkada et ah. 1998) 
















94.00 






MSU 1996 fSwets and Wensr. 1996) 
















-88.50 


-3.00 




Bayesian MAP (Teixeira, 2003) 


1.92 


2.51 


4.18 


6.52 


13.8 


70.15 


67.11 


81.92 






RRGM Standard fRolme. 2003) 


0.59 


0.92 


1.34 


2.42 


9.2 


37.54 


34.26 


88.37 






RRGM Optimised fRolme. 2003) 
















89.80 






PCA MahCosine fCSU. 2003) 


0.84 


1.17 


2.26 


4.43 


10.28 


60.45 


48.90 


85.27 






Gabor features fKeoenekci et ah. 2002) 
















96.30 






Haar+AdaRoost (Jones and Viola, 2003) 








-0.42 


-1.17 






-94.00 


-1.00 




Gabor+AdaRoost (Yane et ah. 2004) 
















-95.20 






SOM (Tan et. ah. 2005) 
















-91.00 






Our hog-Gabor PCA, 900 PCA feat., 3 pt. 


0.08 


0.08 


0.08 


0.08 


0.25 


1.84 


8.92 


98.16 


0.33 


0.33 


Our hog-Gabor PCA, 900 PCA feat., 2 pt. 


0.08 


0.08 


0.08 


0.17 


0.33 


68.81 


24.02 


97.99 


0.33 


15.78 


Our trad. PCA, 900 PCA feat., 3 pt. 


2.76 


4.35 


7.02 


12.63 


38.13 


92.89 


123.40 


82.43 


3.93 


110.99 


Our trad. PCA, 900 PCA feat., 2 pt. 


6.94 


8.61 


13.21 


16.64 


31.44 


99.58 


149.75 


76.90 


5.27 


136.24 


Our hog-Gabor PCA 4x4, 900 PCA feat., 3 pt. 


0.08 


0.08 


0.08 


0.08 


0.08 


0.17 


0.59 


98.49 


0.17 


0.14 



recognition were used 300 featur es. USC 1997 (University of Southern Ca lifor- 
nia) method was developed by (|Wiskott et all 119971 lOkada et all Il998h . For 
recognition they used Gabor Jets and Elastic Bunch Graph Maching (EBGM). 
Faces were resized to 128x128 pixels and normalized using histogram equal- 
ization For recognition were used about 1920 features that correspond to 40 
Gabor filters (5 scales and 8 orientations) at 48 gr aph nodes. MSU 1996 (M ichi- 
gan State University) method was dev eloped bv ISwets and Wene (1996). For 
recog nition they used PCA and LDA. iGSUl (|2003f) Bayes i an MA P (|Teixeira . 
2003| i. EBGM Standard, and EBGM Optimised (|Bolmel . liool face recogni- 



tion methods were developed by the researchers at Colorado State University. 
These methods are similar t o the corresponding methods developed at MIT 
and USC. The ICSUl (|2003h PCA MahCosine method is a traditional PCA 



with c osine -based distance measure between "whitened" feature vectors. ICSU 
(2003) 'or recognition used 130x150 images, faces were masked using ellipti- 
cal mask, unmasked image part was normalized using histogram equalization. 
CSU EBGM method for recognition extracts more than 60 00 features (80 
Gabor features x 80 graph points). iKepenekci et al. (2002) for recognition 
used magnitudes of Gabo r filters and similar distance measures as were used 
by (jWiskott et all 119971 ). For recognition were extracted 40 Gabor features 



(5 scales and 8 orientat ions) and 2 coordinate s of these features at varying 
number of face points. Jones and Viola ( 2003| ) used Haar -like features and 
AdaBoost training. For recognition were used 45x36 images without masking. 
The use small images may be related with the fact that A daBoost training re- 
quires huge computational resources. lYang et al] (2004) used Gabor features 
and AdaBoost training -based recognition method iTan et. al. ( 2005| ) for face 
recognition used Self- Organizing Map (SOM) and soft k nearest neighbor (soft 
k-NN) ensemble method. 



As we can see from the Table 2, the highest first one recognition accuracy 
was achieved by the following methods: our Log-Gabor PCA (98.16% us- 
ing 2-poin t normalization and 98 .49% using 3-point normalization), Gabor 
features (IKepenekci et all l2002h -based method (96.30%), and UMD 1997 
f)Zhao et al!ll998l ) PCA+LDA -based method (96.23%). The bes t EER results 
were achieved by our Log-Gabor PCA (0. 33%), Haar+Ada Boost ([Jones and Viola . 
200l method (-1.00), and UMD 1997 (jZhao et all Il998h PCA+LDA -based 
method (1.09%). It is interesting to note that our traditional PCA with 2- 
poin t normaliza tion achieves lower recognition accuracy than traditional PCA 
of (|CSlAl2003h . But when we combine our traditional PCA with Log-Gabor 
features, our method achieves higher recognition accuracy that many other 
methods. Also we can note that PCA with grayscale features is much more 
sensitive to the chosen image normalization method (76.90% first one recog- 
nition using 2-point normalization and 82.43% using 3-point normalization) 
than our Log-Gabor PCA (98.16% using 2-point normalization and 98.49% 
using 3-point normalization). Using 3-point normalization faces are masked 
and cropped more accurately than using 2-point normalization, and recogni- 
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- O O O O G o o e^e 
75%. rank = 0^2% (11 images) 



o o . 

rec. = 99.58%, rank = 0.59% (7 images) 
rec. = 99.50%, rank = 0.33% (4 images) 



.83%, rank = 0.25% (3 images) 



rec. = 98.41%, rank = 0.17% (2 images) 



rec. = 97.99%, rank = 0.08% (1 image) 



0.08 0.25 0.41 0.58 0.75 0.92 1.08 1.25 1.42 1.59 1.75 1.92 
Rank, 0-100% (of 1196 images; 0.08% - one image) 



FAR-0.01 % FRR-? 85% 



|FAR=0.013%,FRR=2.09% 
J,FAR=0.015%,FRR=1.92% 

FAR=0.02%,FRR=1.76% 

FAR=0.D3%,FRR=1.51% 




FAR=0.05%,FRR=1.26% 

FAR=D07%.FRR=1.00% 
FAR=0.10%,FRR=0 84% 
FAR=0.T1%,FRR=0.67% 
FAR=0.15%,FRR=0.50% 



FAR-FRR-EER-0.33% 
FRR^0.25%,FAR=0 44% 

FRR=025% 

frr=o:di 



FAR=51 .38% 
J%,FAR=B3.23% 



10 



Fig. 6. CMC characteristic of Log-Gabor Fig. 7. ROC characteristic of Log-Gabor 
PCA method PCA method 

tion results using Log-Gabor PCA and 3-point normalization are also better. 
These differences are especially visible when we compare CumlOO, CMCA and 
ROCA values of 2-point ant 3-point methods. It is interesting to note, that 
EBGM -based methods perform positioninig of graph nodes around the face 
also enough accurately and this may be one of the reasons why EBGM -based 
methods and Log-Gabor PCA method using 3-point normalization achieve 
better CumlOO results than other methods that use 2-point normalization. 
Using our method with 3-point normalization in order to achieve 100% cu- 
mulative recognition rate we need to extract from the database only 1.84% of 
images (that is 1196*1.84/100 = 22 images), and using CSU EBGM method 
we need to extract 37.54% of images (that is 449 images). In the last line 
of the Table 2 we present the results of our method if we use 3-point nor- 
malization, 4x4 sliding window (without overlapping), masking, and 19704 
Log-Gabor magnitude features (the number of PCA features remains 900). 
As we can see from these results, using larger number of Log-Gabor features 
(smaller sliding window) we can achieve even better cumulative recognition 
results than using 8x8 window with 2 pixels overlapping. That is in order to 
achieve 99% cumulative recognition rate we need to extract from the database 
only 2 images (0.17%), and in order to achieve 100% cumulative recognition 
rate we need to extract 7 images (0.59%). 

Because in real life situations face recognition methods are usually used with 
automatically detected faces and facial features. So it is desirable to know 
how detection errors affect recognition accuracy, what recognitiom method is 
less sensitive to feature detection errors. Using this information we can de- 
cide how accurately faces and facial features should be detected in order to 
achieve desirable recognition accuracy. In the Table 2 we presented some re- 
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suits of other researchers that used automatical detectio n of faces and facia l 
features (notation "auto"). USC 1997 automatic method ( Okada et all Il998h 
for loacation of face and facial features (eyes, nose, lips, face contour) used 
EBGM with small number of graph nodes (16 nodes). As it is stated by the 
authors, their method locates facial features very accurately, so the difference 
between recognition result s using manually and au tomatically detected fea- 
tures is ~1%. MIT 1996 ( Moghaddam et al. . 1996| ) automatical method for 
detection of eyes and lips used PCA-based detector and probabilistic verifica- 
tion of detected feature locations. Automatical method achieved ~7% lower 
recognition accuracy than the same method that used manual feature detec- 
tion. The results showed that fully automatical USC 1997 method can achieve 
much higher recognition accuracy than MIT 1996 method. But because the 
authors used different methods for detecting faces and facial features and did 
not present any quantitative information about feature detection accuracy, 
we cannot say for sure what recognition method (that was tested without 
automatical detection) it is better to use with automatical feature detection 
method. It is possible that one recognition method is less sensitive to feature 
detection errors than another, but also it is possible that the main differences 
are only in feature detection methods and their feature detection accuracy. 



In order to find out how sensitive is our face recognitiom method to feature 
detection errors and not to bind to concrete feature detection method we man- 
ually shifted the markers of eye centres using pre-defined shift directions and 
distances. For this experiment we used all facial images from the FERET fa 
and fb sets. We shifted only eye markers of fb images using 4 shift directions 
(0, vr/2, 7T, 3vr/2) and the following shift distances: 0%, 2%, 4%, 6%, 8%, 10%, 
12%. Recognition accuracy using each shift distance was calculated as an aver- 
age of 4 results that correspond to 4 directions. Shift distance is calculated as 
a percentage of the distance between manually selected eye centres. For exam- 
ple, if the distance between manually selected eyes is 100 pixels and we wish to 
use 4% shift, then these 4% will correspond to 4 pixels. The results that were 
achieved using 0% shift (it means that no shift is performed and we simply 
use manually selected eye coordinates) were used as a baseline for comparison 
with the results that were achieved using other shift distances. The results 
are presented in Table 3, where notations Firstld and EERd mean absolute 
differences between achieved recognition result using some shift and the result 
without any shift (baseline). For experiments was used 2-point normalization 
method. 



The results (Table 3) showed that our Log-Gabor PCA is less sensitive to fea- 
ture detection errors than traditional PCA and can achieve 89-90% recognition 
accuracy even if one eye is shifted by 10%. In order to create fully automat- 
ical face recognition method and achieve s imilar recognition accuracy that 
was achieved by USC 1997 ( Okada et al. . 1998i) fully automatical method 
we should combine our recognition method with automatical eye detection 
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Table 3 

Face recognition accuracy using shifted markers of eye centres (simulated eye loca- 
tion errors). 



Shift 


Method 




Shifted left 


eye marker 


Shifted right 


eye marker 


size, % 


name 




First 1 


Firstld 


EER 


EERd 


First 1 


Firstld 


EER 


EERd 


0% 


PCA 




76.90 




5.27 




76.90 




5.27 






Log-Gabor 


PCA 


97.99 




0.33 




97.99 




0.33 




2% 


PCA 




73.10 


3.80 


6.05 


0.78 


72.95 


3.95 


6.05 


0.78 




Log-Gabor 


PCA 


97.45 


0.54 


0.40 


0.07 


97.41 


0.58 


0.40 


0.07 


4% 


PCA 




62.45 


14.45 


8.39 


3.12 


61.86 


15.04 


8.56 


3.29 




Log-Gabor 


PCA 


96.76 


1.23 


0.44 


0.11 


96.88 


1.11 


0.50 


0.17 


6% 


PCA 




47.78 


29.12 


12.32 


7.05 


45.04 


31.86 


12.59 


7.32 




Log-Gabor 


PCA 


95.65 


2.34 


0.65 


0.32 


95.31 


2.68 


0.63 


0.30 


8% 


PCA 




31.17 


45.73 


16.38 


11.11 


31.30 


45.60 


17.22 


11.95 




Log-Gabor 


PCA 


93.41 


4.58 


0.96 


0.63 


92.95 


5.04 


0.96 


0.63 


10% 


PCA 




19.21 


57.69 


21.51 


16.24 


19.29 


57.61 


22.53 


17.26 




Log-Gabor 


PCA 


90.02 


7.97 


1.46 


1.13 


89.02 


8.97 


1.55 


1.22 


12% 


PCA 




11.40 


65.50 


26.46 


21.19 


11.92 


64.98 


27.13 


21.86 




Log-Gabor 


PCA 


84.29 


13.70 


2.24 


1.91 


83.24 


14.75 


2.36 


2.03 



method that detects the centres of eyes with smaller than 6% shifts (total 
shift for both images) when compared to manually selected eye centres. Those 
readers who are interested in automatical face and f acial features dete ction 
methods can find an overview of such methods in ((Yang et all l2002h and 
flPerlibakasl boOSl b 



In real life situations the accuracy of face recognition is also affected by many 
other factors like aging and manual change of appearance (hairstyle, makeup), 
image input (camera position) and illumination conditions. It is natural that 
we cannot have the same looking faces and the same imaging conditions after a 
longer time period. In order to find out how recognition accuracy is affected by 
these factors we performed recognition experiments using the following probe 
sets of images from the FERET database: dupl - 722 images of 243 persons, 
at least 2 images with different expressions per person, time interval from 
fa images is 0-34 months, photographs of 166 persons were taken after some 
period of time (not the same day than fa images); dup2 - 234 images of 75 
persons, at least 2 images with different expressions per person, time interval 
from fa images is more than 18 months; and fc - 194 images of 194 persons 
that were acquired on the same day, but with different camera position and 
illumination. 



The results (Table 4) showed that our Log-Gabor PCA method achieves 
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Table 4 

The influence of aging and illumination to face recognition accuracy using 1196 fa 
gallery images (FERET database) and the following probe sets: dupl, dup2, fc. 



Method and its authors 


First one rec, 


0, 100%] 


EER, [0, 100%] 




dupl 


du P 2 


fc 


dupl 


dup2 


fc 


MIT 1996 (NIST, 2001) 


57.60 


34.20 


32.00 


17.70 


21.20 


18.00 


MIT 1996 auto fPhillios et al.. b. 19981 


~50.00 


- 


- 


- 


- 


- 


UMD 1996 (ThilliDs et al.. b. 19981 


~32.00 


~9.00 


-30.00 


- 


- 


- 


UMD 1997 CNIST. 2001) 


47.20 


20.90 


58.80 


12.60 


13.40 


10.00 


MSU 1996 fPhillins et al.. b. 19981 


~33.00 


~17.00 


~32.00 








USC 1997 (NIST, 2001) 


59.10 


52.10 


82.00 


13.30 


14.20 


5.10 


USC 1997 COkada et al..J^98l 


62.00 


52.00 


82.00 








USC 1997 auto fOkada et al.. 1998) 


61.00 


52.00 


80.00 








Gabor feat. (KeDenekci et al.. 2002) 


58.30 


47.40 


69.60 








Our trad. PCA, 900 PCA feat., 2 pt. 


44.74 


35.04 


62.89 


13.99 


19.03 


9.79 


Our Log-Gabor PCA, 900 PCA feat., 2 pt. 


72.44 


65.81 


90.21 


3.60 


4.70 


1.03 



8-10% higher recognition accuracy and at least 4% lower EER than other 
compared methods. Our recognition results showed that even using single 
training image per person we can improve recognition accuracy of face images 
that were took after longer time period. But also it is obvious that different 
imaging conditions after longer time period significanly reduce face recognition 
accuracy of all compared methods and that for such difficult tasks we need 
better image normalization and feature extraction techniques. 



We also perform ed several face recognition experiments using the AR database 
f Martinea . 1998). This database was created by A. Martinez and R. Benavente 
at Computer Vision Center, Purdue University in 1998. It contains facial pho- 
tographs of 126 persons with strictly controlled facial expressions and lighting. 
The size of images is 768x576 pixels. Images of each person were captured 
in two sessions (si, s2) that were separated by two weeks time. From this 
database we used 1610 images of 115 persons (14 images per person = 2 ses- 
sions x 7 images per session). We used the following images: neutral (ne), 
happy (ha), angry (an), and screaming (sc) expressions; neutral expression 
with left illumination source (lis) turned on, right illumination source (ris) 
turned on, and both illumination sources (bis) turned on. For training and 
as a galery set we used 115 images with neutral expression from the first 
session (sine). For recognition we used the following probe sets that corre- 
spond to different type of transformation (neutral, expression, illumination) 
f)Wang and TansJ . |2003) : slexpr (slha, slan, slsc images), slillum (sllis, slris, 
slbis images), s2neutral (s2ne images), s2expr (s2ha, s2an, s2sc images), and 
s2illum (s21is, s2ris, s2bis images). First one recognition results using these sets 
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are presented in Table 5. The last lines of this table also present EER results 
of our methods. For experiments was used 2-point normalization method. 



Table 5 

Face recognition results using AR database. 



Method and its authors 


Firs one 


rec. results using different probe sets 




slexpr 


slillum s2neutral 


s2expr 


s2illum 


r \ vvanE ana lang, zuuoj 
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04.4 


7 

DO. 1 
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HiDKjrLvi ijaDor iearures ( vvanE ana laiiE, zuuoj 




oD. 1 


P.P. 7 
DD. i 


9 


T^Rf^A/T PoVtar foolnroc _1_ 

xijo^jivi v^raDor ieatures -t- 










odyub iiiciiciiiiig, \ vvdiiii diiu idiiti, zuuoj 




ao.o 


oD.U 


Rfi 7 
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PCA + optical flow (Martinez. 2Q^^) 


83.00 








Motion estimation (Martinez, 2003. b) 


84.67 








Our traditional PCA, 100 feat., 2 pt. 


70.43 


62.90 92.17 


58.52 


46.67 


Our Log-Gabor PCA, 100 PCA feat., 2 pt. 


85.51 


82.90 99.13 


77.39 


63.48 






EER results of our 


methods 




Our traditional PCA, 100 feat., 2 pt. 


6.96 


5.51 2.61 


11.02 


10.72 


Our Log-Gabor PCA, 100 PCA feat., 2 pt. 


3.48 


3.19 0.87 


6.67 


6.67 



The results (Table 5) showed that our method achieves slightly higher recog- 
nition accuracy than other compared methods that used single training image 
(with neutral expression) per person. It is important to note that the compar- 
ison of different methods in the Table 5 is not very exact, b ecause for exper- 
imen ts different authors used diff erent number of i m ages: (IWang and Tane . 
l2003h used images of 90 persons, (jMartir1e7l . Enn3.ah . (jMartiue7l . l2nn3.bh used 
images of inn persons, and we used images of 115 persons. The results showed 
that EBGM Gabor features ( features are extracted a t graph nodes) and Bayes 
matching -based algorithm (jWang and Tane , 2nn3h can achieve much higher 
recognition accuracy than our method and all other methods (our method 
achieved better results only when recognizing faces with neutral expressions). 
But for training of this method we need multiple images per person, and this 
method was trained using 7 images per person from the first session with differ- 
ent expressions and illumination conditions. All other compared methods for 
training used single image per person (image with neutral e xpression f rom the 
first session). It is interes ting to no t e PCA + optical flow (jMartinezl . bnn3.ah 
and Motion estimation ( Martinet . l2003.tf ) -based recognition methods were 
specially designed for recognizing expression -variant faces, and these methods 
use weigting of facial features in order to reduce the influence of changed ex- 
pression to the accuracy of face recognition. Perhaps these weigting methods 
could improve recognition accuracy of our face recognition method, and in 
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the future we are going to investigate different feature weighting and masking 
methods in order to improve recognition accuracy of expression -variant faces. 



7 Conclusions and future work 

In this article we proposed a novel face recognition method based on Principal 
Component Analysis (PCA) and Log-Gabor filters. The experiments showed 
that using the proposed combination of Log-Gabor features and sliding win- 
dow -based feature selection method, Principal Component Analysis, "whiten- 
ing" , and cosine -based distance measure we can achieve very high recognition 
accuracy (97-98%) and low error rates (0.3-0.4% Equal Error Rate) using the 
FERET database that contains photographs of more than 1000 persons. The 
results of our algorithm are among the best results that were ever achieved 
using this database. In the future we are going to investigate the possibilities 
of using decomposed Log-Gabor feature vectors and multiple PCA spaces in 
order to have the possibility of using this method with an unlimited number 
of training images. Because the results of all compared methods showed that 
the accuracy of face recognition is very affected by the lighting conditions, in 
the future we are going to investigate different lighting normalization methods 
and test them with the Log-Gabor PCA face recognition method. 
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